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PJ ■ ABSTRACT 

^ ' It is shown that the cuspy density distributions observed in the cores of ehiptical 

galaxies can be reahzed by dissipationless gravitational collapse. The initial models 
fSJ ' consist of power-law density spheres such bs. p (xr~^ with anisotropic velocity disper- 

Osl I sions. Collapse simulations are carried out by integrating the collisionless Boltzmann 

equation directly, on the assumption of spherical symmetry. From the results obtained, 
the extent of constant density cores, formed through violent relaxation, decreases as 
^ \ the velocity anisotropy increases radially, and practically disappears for extremely ra- 

dially anisotropic models. As a result, the relaxed density distributions become more 
cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity 
anisotropy could be a key ingredient for the formation of density cusps in a dissipa- 
(^ ■ tionless collapse picture. The velocity dispersions increase with radius in the cores 

Qv^ I according to the nearly power-law density distributions. The power-law index, n, of 

^■^ . the density profiles, defined as p oc r~", changes from n « 2.1 at intermediate radii, to 

a shallower power than n w 2.1 toward the centre. This density bend can be explained 
from our postulated local phase-space constraint that the phase-space density acces- 
sible to the relaxed state is determined at each radius by the maximum phase-space 
O ' density of the initial state. 

H ; 

"^ \ Key words: galaxies: formation - galaxies: kinematics and dynamics - galaxies: 

C^ . structure - methods: numerical. 
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H 1 INTRODUCTION (1997) have pointed out from their high resolution simu- 

. - - ■ rr 7 , 7 rr /^ latlon wlth about 8 x 10^ particles that the resulting cen- 

Recent Hubble Space Telescope observations (Crane et al. . , , •. m r i i .. i i . .i 

, „„„ ^ 7 , „„ , T rr- 7 1 „„ - r 7 ^ral deusity pronles oi dark matter haloes are steeper than 

1993; lerrarese et al. 1994; Jane et al. 1994; Lauer et al. _i . , . 7.77.. . 7,7, 

.„ , „ ^ ,7 7 7 ' ^7 7 pocr m a hierarchical clustermg scenario, and that merg- 

1995; Byun et al. 1996; Gebhardt et al. 1996; Faber et al. f^ , ■^, .-rr ^ I- j- ■ 

^^^rj\ 1, 1 J iu 4- ir_4-- 1 1 • J 1 4- T, „ ^'^S processes among clumps with diilerent binding energies 

would play a significant role in the formation of such steep 

density cusps. Moore et al. (1998) have also demonstrated 

that very high resolution simulations with nearly 3 x 10® 

particles yield steep inner density profiles with p oc r"^* 

(see however Kravtsov et al. 1998). On the other hand, 

p (X r with 0.5 ;$, ?^ ^ 1 tor bright elliptical galaxies, , ,, 7 .^ 7 7 7 r 7 ^7 7 

'^ 7 7 7 7 shallower density cusps could be tormed through merging 

and the other is a steep cuspy core described by approxi- , , , . 7 . 7 . 7. . 7 77 7 • ,7 

_2 ^ ^ . ,, J . o T^ . 7 , „„^N between galaxies which individually have a massive central 

mately p oc r tor tamt ones JVlerritt & Iridman 1996 . , , , , , , 7 m- 7 • T-n ■ 7 • /ir.r>^\ n ^7 • 

„ , , , , , , 7777 black hole, as shown by Makmo & Ebisuzaki (1996). For this 

Several cosmological simulations, based on a standard cold , ■ j, r ■ tvt? otvti- f-,nnn\ 7 7 • j 

° kind ol merging, Nakano & Makmo (1999) have explained 

dark matter scenario, have demonstrated that cuspy density , , ,7,7 77 7 .777777 

_, , why a galaxy that has swallowed a massive black hole leaves 

profiles with n oc r might be a natural end-product or hi- 7 ,7 . ,7 , tt .77.7 7 

/T^ 7 7 o ^ 77 TVT Bi shallow cusp m the centre. However, the physical mecha- 

erarchical clustering (Dubmski & Carlberg 1991; Navarro, . , , . ,.,, . , , 

„ ,„,,„ ,„ ,„,, „,,, nism to generate steep cusps still remains unsolved. 

Frenk & White 1996, 1997). Recently, Fukushige & Makino ^ f f 

From the simulations mentioned above, we might con- 
clude that clumps are needed to produce steep cusps while 
* E-mail: hozumi@sue.shiga-u.ac.jp (SH); merging between galaxies with massive central black holes 

burkert@mpia-hd.mpg.de (AB); fujiwara@kcua.ac.jp (TF) is required to generate shallow cusps. However, it might be 
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1997) have revealed that elliptical galaxies do not have con- 
stant density cores but have cusps which continue toward 
the centre until the resolution limit. In addition, it has 
been found that the density profiles of the cores are divided 
into two types: one is a shallow cuspy core represented by 



2 S. Hozumi, A. Burkert and T. Fujiwara 



too early to abandon a simple dissipationless collapse pic- 
ture, which includes neither dumpiness nor massive central 
black holes, for the formation of cusps. We know that the 
collapse of uniform-density spheres results in large constant 
density cores (van Albada 1982; Fujiwara 1983), while lit- 
tle is known about the details of the collapse of power-law 
density spheres. Of course, there are cases for the collapse 
of non-uniform density spheres. For example, Hozumi, Fuji- 
wara & Kan-ya (1996) have carried out collapse simulations 
of power-law density spheres only in the initial contraction 
phases in order to demonstrate the predominance of the 
growth of the tangential velocity dispersion over that of the 
radial velocity dispersion until the most contracting phase. 
In addition, Londrillo, Messina & Stiavelli (1991) examined 
the collapse of generalized Plummer models; their models 
have large cores, and steep density haloes with p oc r~^ or 
p oc r"'^, so that the collapse behaviour over most parts of 
the system is considered similar to that of a uniform-density 
sphere. Burkert (1990) paid attention to degenerate cores 
arising from the dissipationless collapse of spherical stellar 
systems with different initial density gradients and virial ra- 
tios. In another instance, Cannizzo & HoUister (1992) stud- 
ied the relation between the power-law index of initial den- 
sity distributions and the final axis ratios of the systems. 
However, the structures at sufficiently small radii where the 
density cusps are observed in real galaxies were not a main 
concern in these examples. 

Although many dissipationless collapse simulations 
have been carried out so far, there has been no sufficiently 
high resolution to determine fine structures at very small 
radii, such as density cusps. In addition, the systems in 
numerical studies were assumed to have isotropic velocity 
distributions at the beginning. That is, the initial models 
for dissipationless collapse simulations have never included 
anisotropic velocity distributions, despite the fact that stel- 
lar systems can support velocity anisotropy unlike gaseous 
systems. 

In this paper, we examine the influence of initial ve- 
locity anisotropy on the relaxed density profiles in the core, 
assigning enough resolution to analyse the innermost region, 
and show that radially anisotropic models could produce dif- 
ferent cuspy density profiles. In Section 2, we describe the 
initial models in which anisotropic velocity distributions are 
taken into account, and the numerical method. In Section 3, 
the results of collapse simulations are presented. In Section 
4, we explain the resulting density profiles from a viewpoint 
that the maximum phase-space density accessible to the re- 
laxed system is determined at each radius by that of the 
initial system. Conclusions are given in Section 5. 



2 MODELS AND METHOD 

We carry out collapse simulations of power-law density 
spheres on the assumption of spherical symmetry. The initial 
virial ratios, -q = 2T/\W\, used are rj = 0.1 and -q — 10~^'^, 
where T is the kinetic energy, and W is the potential energy. 
According to existing simulations (Polyachenko 1981, 1992; 
Merritt & Aguilar 1985; Min & Choi 1989; Aguilar & Merritt 
1990; Londrillo et al. 1991; Cannizzo & HoUister 1992; Udry 
1993) like those studied here, our models would be affected 
by the radial orbit instability if spherical symmetry were not 



assumed, because of the small initial virial ratios. However, 
Aguilar & Merritt (1990) have shown that the origin of the 
radial orbit instability is different from that of the universal 
density profile observed in elliptical galaxies, that is, the de 
Vaucouleurs R^'"^ law in projection (see also Burkert 1993). 
Therefore, we are allowed to assume spherical symmetry as 
long as we focus only on the spherically averaged density 
profiles and related physical quantities of collapsed objects. 
The velocity anisotropy is introduced into the initial 
models by adopting the distribution function, /, as 
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where po is the density, Ur is the radial velocity dispersion, 
and (Tx is the tangential velocity dispersion with r, u, and 
j being the radius, radial velocity, and angular momentum, 
respectively. We begin with a moderate power-law density 
profile given by 



po(r) = {M/2nRl)r- 



(2) 



where M is the total mass and Rq is the radius of the 
sphere. We measure the velocity anisotropy by a = 2Tr/T±, 
where Tr and T± are the kinetic energies in radial and tan- 
gential motions, respectively. In our models, a is equiva- 
lent to 2(jJ;/(Jx, as is found from equation (|l|). We take 
a — 0.2,0.5,1,2,6, and 10 for each value of the virial ra- 
tios. 

We use Fujiwara's (1983) phase-space solver which in- 
tegrates the collisionless Boltzmann equation directly, with 
the help of a splitting scheme (Cheng & Knorr 1976) . Since 
this method leads to the results that are not subject to ran- 
dom fluctuations, we can obtain smooth density and velocity 
dispersion profiles. We employ {Nr, Nu, Nj) = (300, 301, 50), 
where Nr, Nu, and A'j are the numbers of mesh points along 
the radius, radial velocity, and angular momentum, respec- 
tively. A logarithmic grid is used in the radial direction. An 
equal interval grid is assigned to the radial velocity while 
a quadratic interval is adopted for the angular momentum. 
The units of mass, M, and the gravitational constant, G, 
are taken to be Af = G = 1. The unit of length is deter- 
mined from the relation such that rjx Ro — 1. This choice of 
Ro leads to approximately the same core sizes for uniform- 
density spheres with isotropic velocity distributions (Fuji- 
wara 1983). As a characteristic time-scale for the collapse, 
we adopt the free-fall time, defined by tff = nWR^/SGM. 
We obtain tg = 35.1 for 77 = 0.1 and tft = 197 for tj = 10"^'^ 
in our system of units. The minimum and maximum radii 
assigned to grid points are set to be -Rmin ~ 0.01 and 
J^max = 40.0, respectively, for 77 — 0.1, and -Rmin = 0.01 and 
-Rmax = 49.9, respectively, for 77 = 10~^'^. Test simulations 
have shown that resulting density and velocity dispersion 
profiles at small radii do not depend on the replacement of 
Rmin ~ 0.01 by -Rmin ~ 0.005 with morc finely divided mesh 
points in the angular momentum. 

Hozumi & Hernquist (1995) and Hozumi (1997) traced 
the orbits of stars on grid points backward to t = ev- 
ery suitable time step in order to avoid the numerical diffu- 
sion generated by the repeated interpolation required by the 
splitting scheme. However, we do not adopt this method in 
our simulations. This is because we are interested not in the 
details of the phase-space distributions but mainly in low- 
order moments of the distribution function such as density 
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and velocity dispersion, in addition to saving computational 
time. 



bution. In reality, as shown by Fig. 2, the imperfect linearity 
of the velocity dispersion against radius in logarithmic scale 
does not realize a perfect power-law density distribution. 



3 RESULTS 

We stopped the simulations at f = 120 for 77 = 0.1 and 
at £ = 500 for 77 — 10~^'^. These times correspond to about 
3 iff. They are, nevertheless, sufHciently long that the result- 
ing cores have relaxed completely. Since our main concern 
is the core regions, these times satisfy our requirements for 
analysing the relaxed density and velocity dispersion pro- 
files near those regions. There were practically no escapers 
for rj = 0.1, while about 20 per cent of the total mass with 
positive energies expanded beyond -Rmax and escaped for 



V = 10" 



Consequently, the total energy was conserved 



to better than 0.39 per cent for all the models with rj = 0.1 
except for the model with a = 0.2 in which the total energy 
within 7?max changed by 1.8 per cent owning to a small frac- 
tion of escapers. On the other hand, the total energy within 
-Rmax changed by at worst 10.7 per cent for all the models 
with 1] = 10~^'^ because of a large fraction of escapers. 

In Figs, la and lb, we show the resulting density pro- 
files for 7] = 0.1 and 77 = 10~^'^, respectively. The half-mass 
radii, rh, for rj — 0.1 are about 2.6 to 3.2, while they are 
about 7.8 to 9.7 for 77 — 10~^'^. Thus, our adopted mini- 
mum radii correspond to -Rmin ~ 0.003 rh for 77 = 0.1 and 
^min ~ 0.001 Th for T] = lO"^"'. If the half-mass radius is 
regarded practically as the effective radius, r^, of the de Vau- 
couleurs law, the minimum radii used are sufficiently small 
and correspond to roughly 10 pc, because re is typically 3 
kpc for bright ellipticals and is smaller for faint ones (Kor- 
mendy 1977). We point out that most galaxies are not well 
resolved inside to 10 pc according to the data cited in the 
paper of Gebhardt et al. (1996). 

We can see from Fig. 1 that the density distributions 
at small radii become more cuspy as a increases, and that 
their slope becomes steeper with increasing a for both virial 
ratios. As a result, the core region, specified by a nearly 
constant density distribution, decreases with increasing a. 
In particular, it is to be noticed that there is no constant 
density core for the most radially anisotropic models, that 
is, the models with a — 10. On the other hand, the den- 
sity profiles at intermediate radii are well-approximated by 
p oc r~^'^ , regardless of the virial ratios and the velocity 
anisotropy. 

The relaxed velocity dispersion profiles for 77 = 0.1 and 
T] = lO""*^'^ are shown in Figs. 2a and 2b, respectively. These 
figures indicate that the relaxed velocity distributions no 
longer remain anisotropic in the inner region of the cores 
even though we begin with a model which has a large ve- 
locity anisotropy at all radii. Thus, it seems very difficult to 
produce and maintain velocity anisotropy in the cores only 
by gravitational collapse. The isotropic region that roughly 
corresponds to the core becomes smaller, as the velocity 
anisotropy increases. It is found from Figs. 2a and 2b that 
the density of the inner regions necessarily decreases with ra- 
dius because the velocity dispersions rise with radius in the 
cores. In particular, for a = 10, the velocity dispersions in 
the cores increase almost linearly with radius in logarithmic 
scale. Then, the hydrostatic equation requires that the cor- 
responding density profiles should follow a power-law distri- 



4 DISCUSSION 

4.1 Phase-space constraint 

CoUisionless stellar systems suffer the constraint that the 
phase-space density accessible to the relaxed state is less 
than the maximum phase-space density of the initial state. 
On the basis of this phase-space constraint, we can explain 
the density bend which occurs around the edge of the cores. 

As far as the core generated through violent relaxation 
is concerned, the phase-space density does not decrease sub- 
stantially from the initial value at least for spherically sym- 
metric systems (Fujiwara 1983; Hozumi 1997). Our adopted 
initial virial ratios are small enough to induce a radial orbit 
instability, unless spherical symmetry is imposed. Neverthe- 
less, it may be unlikely that the phase-space density in the 
core suffers a substantial decrease after violent relaxation 
according to the results of May & van Albada (1984): they 
have shown that the phase-space density in the core is al- 
most conserved after violent relaxation even for aspherical 
collapses, although they examined this behaviour only for 
initially homogeneous density distribution models. The ef- 
fects of aspherical collapses through the radial orbit insta- 
bility are left to be investigated as a future problem. 

Here, admitting the non-decreasing nature of the phase- 
space density in the core, we postulate a local phase-space 
constraint that the maximum phase-space density in the re- 
laxed system is determined by the initial conditions at each 
radius in the core regions for power-law density spheres. 
From equation (hi), the local maximum phase-space density 
at the beginning, /m, is defined by 



Mr) 



Po{r) 



^7r3/2cr,cr2 



(3) 



If the initial density distribution is given by a power law 
like po oc r~", the local maximum phase-space density of 
the initial system is represented in terms of the anisotropic 
parameter, a, and the virial ratio, rj, by 



fm{r) = 



T-l/2 



— n —n 

T} r 



(4) 



8\/27r5/2 (3-71)1/2 

Notice that equation (W) holds even though a depends on 
radius. The cumulative mass within r is M{r) = (r/i?o)'^~"'- 
For the special case where a is independent of r, we find the 
cumulative mass with phase-space density greater than /m 
by solving equation Wm with respect to r, 



M{> U) 



(5 - 27l)3(3-n)/(2n) 



[8^7^5/2] (3-n)/n (3 _ ^) {3-n)/(2n) 
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This equation shows that M(> /m) is independent of rj, 
because we have determined Ro from the relation such that 
77 X _Ro = 1. Thus, our discussion described below holds 
whatever value the virial ratio will be. From equation (5), 
we obtain the relation between M and fm as 
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Since our initial models have po cc r~ , equation (|6|) pre- 
dicts the relation M{> /m) oc /m^. The resulting M{> fm) 
after violent relaxation cannot exceed M{> fm) of the 
initial system. As the relaxed density distributions show 
p oc r~ '^ at intermediate radii, regardless of the values of a, 
M{> fm) cc fm would hold approximately at such radii. 
As illustrated in Fig. 3, for large values of fm there un- 
avoidably emerges a region in which M{> fm) for p ex r" '^ 
becomes larger than that for p cc r~^ . Since large fm corre- 
spond to small r as indicated by equation (M), we find from 
Fig. 3 that the density profile with p cc r~ '^ cannot con- 
tinue down to sufficiently small radii owning to the lack of 
phase-space density. Therefore, the density bend necessarily 
arises as long as the initial density distribution is shallower 
than Po oc r~^'^ . 

As a piece of evidence that the local phase-space con- 
straint is valid, we present in Fig. 4 the plots of phase par- 
ticles on the radial mesh points at the final state for 77 = 0.1 
with a = 0.2, 1, and 10 in the (/,r)— plane. This figure re- 
veals that phase particles with high phase-space density are 
not scattered in a wide range of radius. Thus, at least in the 
central region, mixing of high phase-space density particles 
with low phase-space ones does not occur efficiently through 
violent relaxation. If we pick out the maximum phase-space 
density at each radius from that distribution of the phase 
particles which is shown in Fig. 4, we can obtain, in a prac- 
tical sense, the local maximum phase-space density, fm, in 
the relaxed state. Fig. 5 is depicted in such a way to show 
the relaxed fm against r near the central region. We can see 
from this figure that the relaxed fm is well-approximated 
by fm oc r~^ at r ^ 0.1, although the slope of the model 
with a = 10 is closer to —1.6 rather than —2. The rela- 



tion fir 



is naturally derived from equation (M) for 



the density profile p oc r~^ which is in fact obtained at in- 
termediate radii in our simulations, though it is precisely 
p oc r~^'^. Therefore, our explanation for the density bend 
based on Fig. 3 is justified qualitatively. 

We now turn to the tendency that larger values of a 
result in more cuspy density distributions. When we be- 
gin with po oc r~^, equation (H) indicates that M(> fm) 
is proportional to (a -|- 2)^/q. This function has a minimum 
at a = 1. Hence, M{> fm) increases with increasing a if 
a > 1. This means that M(> fm) for p cc r~^ intersects 
with M(> fm) for p oc r~^ with a = 10 at a larger fm than 
that with a — 1 (see Fig. 3). Thus, p oc r~^'^ can continue 
down to smaller radii with increasing q, provided that the 
relaxed density profiles at intermediate radii do not change 
greatly among the models with different values of a, as is 
demonstrated by our simulations. In this way, we can under- 
stand that the large a models produce nearly cuspy density 
profiles down to sufficiently small radii. 

On the other hand, M{> fm) decreases with increas- 
ing a if Q < 1. Thus, if a < 1, smaller values of a should 
also, in principle, lead to a density bend closer to the centre 
(see Fig. 3). In reality, this is not the case, as illustrated in 
Fig. 1. This can be understood from Fig. 6 which shows the 



fractional mass, dM{j)/dj, with the angular momentum, j, 
for the initial models. These fractional mass distributions 
remain unchanged throughout the evolution because spher- 
ical symmetry is assumed. Fig. 6 shows that the mass of 
stars with small values of the angular momentum decreases 
as Of (< 1) decreases for j ^ 0.2. If the p (x r~^ profile could 
continue down to the very centre, and if most stars would 
be on circular orbits, dM{j)/dj would be equal to V47rp, , 
where p, is the reference density at r = 1. The values of 
the imagined dM{j)/dj are 0.45 for 77 = 0.1 and 0.26 for 



7? = 10~ 



For these values of dM{j)/dj, Fig. 6 indicates 



that the smaller a models contain less stars with small j. 
Since the stars with small values of the angular momentum 
can pass close to the centre, the deficiency of small j stars 
means that the smaller a models cannot populate a suffi- 
cient number of stars down to smaller radii. Therefore, the 
radius of the density bend moves outward as a decreases 
if a < 1, even though there is sufficient phase-space den- 
sity available. If spherical symmetry were not imposed, some 
degree of mixing in angular momentum could be expected 
(May & van Albada 1984), and so, p oc r~^'^ might continue 
down to smaller radii as a decreases when q < 1. 

We can infer from Figs. 1 and 6 that the models with 
a > 1 would not suffer the limitation, described above, on 
the central density distribution arising from the amount of 
small angular momentum stars. This is because we can ex- 
plain the shift of the density bend toward smaller radii with 
increasing a simply from the local phase-space constraint 
like that illustrated in Fig. 3. 

We have found that the velocity anisotropy could be ef- 
fective in the formation of cuspy density profiles. However, 
it is unclear how the system can acquire such a large radial 
anisotropic velocity dispersion. It may be difficult to real- 
ize a large velocity anisotropy in each stellar system before 
collapse. In a cosmological situation based on cold dark mat- 
ter, clumps are first formed and then collapse to merge into 
a large system. Thus, a system consisting of clumps may 
have anisotropic velocity dispersions. Therefore, it will be 
important to study the velocity distribution just after such 
a system has decoupled from the Hubble flow. 

4.2 Relation between initial and final density 
powers 

We have found that more radially anisotropic velocity dis- 
tributions lead to more cuspy density profiles by gravita- 
tional collapse. Here, we examine how seriously the initial 
power index of density distributions influences the final one. 
We carry out collapse simulations of density distributions 
p oc r~" with n—0.5, 1, and 2, again on the assumption of 
spherical symmetry using the collisionless Boltzmann code. 
The anisotropic parameter and virial ratio are chosen to be 
Q = 10 and Tj — 0.1, respectively. The other numerical pa- 
rameters are the same as those used in Section 2. The final 
density profiles are shown in Fig. 7. The half-mass radii, 
Th, are about 2.3 to 3.2. We can see that the density dis- 
tributions at intermediate radii are well-approximated by 
p oc r~^'^, and that initially steeper density distributions 
result in steeper cuspy density distributions, although the fi- 
nal power-law indices within the radius of the density bends 
become always shallower than the initial ones. In particu- 
lar, the density power for the model with p oc r~^ remains 
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almost unchanged after the collapse, so that practically no 
density bend is found. This fact is in good agreement with 
the results of Burkert (1990), and proves the validity of the 
local phase-space constraint that we have postulated in the 
previous subsection, because the relaxed density distribu- 
tion with p ex r~^'^ can almost continue toward the very 
centre without the phase-space constraint if the initial den- 
sity distribution is p ex r^^, as indicated by Fig. 3. 

Taking into consideration our results obtained here, the 
steep cusps with p ex r^^ observed in faint ellipticals could 
originate from the fluctuation spectrum that realizes density 
distributions with p ex r^^ . According to observations, low- 
luminosity ellipticals that have the steep cusps are isotropic 
in the sense that they show iy/a)* ~ 1, where {V/a)* is 
the ratio of the rotation parameter V/a to the value for an 
isotropic oblate spheroid flattened by rotation, V is the max- 
imum rotation velocity, and g is the mean velocity dispersion 
inside one-half of the effective radius (Kormendy & Bender 
1996; Faber et al. 1997). In addition, these ellipticals are 
rapidly rotating. Unfortunately, we cannot discuss iy/a)* 
values for our end-products, because our models include no 
net rotation. Therefore, we will need to investigate the ef- 
fects of rotation on the steep cusps using three-dimensional 
collapses in order to make a detailed comparison with ob- 
servations. 

The resulting power-law indices of the density distribu- 
tions within the radius of the density bend correlate with 
those of the initial density distributions. The difference in 
power-law index of the cusps originates from a difference in 
that of initial density distributions. In cosmological simu- 
lations with high resolution, based on a standard cold dark 
matter scenario, the final density distributions show a rather 
steep density cusp, such as p oc r^^'* (Fukushige & Makino 
1997; Moore et al. 1998). Probably, this kind of violent relax- 
ation would not produce shallow cusps observed in bright el- 
lipticals, so that we might require a merger between galaxies 
each of which contains a massive central black hole (Makino 
& Ebisuzaki 1996; Nakano & Makino 1999). However, ff 
the circumstances for the formation of bright ellipticals are 
somehow very different from those of faint ellipticals, it is 
conceivable from our results that bright ellipticals could have 
been formed from the fiuctuation spectrum that corresponds 
to p oc r^" with 0.5 ^ n ^ 1 approximately. 



5 CONCLUSIONS 

Collapse simulations of power-law density spheres such as 
p oc r~^ with anisotropic velocity distributions have been 
performed on the assumption of spherical symmetry. We 
have found that the resulting density profiles are well- 
approximated by p (X r~^'^ at intermediate radii, regard- 
less of the velocity anisotropy and the virial ratio. However, 
this universal density profile cannot continue down to suffi- 
ciently small radii, and so, the density bends necessarily. We 
have explained this density bend from the local phase-space 
constraint that the phase-space density accessible to the re- 
laxed system is determined at each radius by the maximum 
phase-space density of the initial system. 

As the radial velocity anisotropy increases, the density 
profiles of the end-products after violent relaxation become 
more cuspy. In particular, the extremely radially anisotropic 



models with q = 10 have practically no fiat cores. Thus, our 
results imply that the velocity anisotropy could play a sub- 
stantial role in the formation of the density cusps observed in 
elliptical galaxies as long as we rely on a simple dissipation- 
less collapse picture. In reality, it will be difficult to produce 
a large velocity anisotropy in individual systems. Concern- 
ing this problem, however, we point out that a system com- 
posed of clumps, as can be envisaged in the standard cold 
dark matter scenario, might have a large velocity anisotropy 
when it has decoupled from the Hubble expansion. In addi- 
tion, we have examined the relation between the initial and 
final power-law indices of the density profiles by performing 
collapse simulations for p ex r~" models with n=0.5, 1, and 
2. In these cases, the resulting density profiles at intermedi- 
ate radii are also well-approximated by p ex r^^'^ , regardless 
of the initial density power indices. The results show that 
iintially steeper density profiles result in steeper cuspy den- 
sity profiles. Therefore, the difference in the power index of 
the cusps between faint and bright elliptical galaxies might 
arise from the difference in the power index of initial density 
profiles. In particular, the steep cusps with p ex r~^ observed 
in faint ellipticals could result from initial density distribu- 
tions with p ex r~^ , because such density distributions can 
avoid the local phase-space constraint. 

As a by-product, we have also found that the resulting 
cores no longer show velocity anisotropy even though we 
start with a model which has a large velocity anisotropy 
at all radii. This implies that simple gravitational collapse 
could have difficulty in generating and maintaining velocity 
anisotropy in the cores. 

Our results shown here are obtained on the assumption 
of spherical symmetry. If aspherical collapses are allowed, 
the radial orbit instability that arises from small initial virial 
ratios might affect the density distributions of end-products. 
The study of three-dimensional collapses is in progress. 
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Figure 1. Relaxed density profiles for the models with (a) 2T/\W\ = 0.1 and (b) 2T/\W\ = 10 ^■^. The straight dotted lines show 
the density profiles with p oc r~^-^. The half-mass radii, rj^, for 2T/\W\ = 0.1 are about 2.6 to 3.2, while they are about 7.8 to 9.7 for 
2T/\W\ = 10-1-5. 
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Figure 2. Relaxed velocity dispersion profiles for the models with (a) 2T/|VK| = 0.1 and (b) 2T/\W\ = 10 ^'^. The solid lines denote 
the square of the radial velocity dispersion, and the dashed lines show that of the tangential velocity dispersion. 
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M(> /m) 




Figure 3. Cumulative mass, M{> /m), with phaso-spaco density greater than /m against the maximum phase-space density of the initial 
models at each radius, /m- The solid lines denote the models with a = 0.2, 1, and 10 for the initial density distributions, p oc r~^. The 
dotted line represents the mass distribution for p ccr~^ . Notice that the mass distribution in phase-space density is independent of the 
initial virial ratio because the initial size of the sphere, Ro, is set to be a reciprocal of the virial ratio. In this plot, large /m correspond 
to small r. 
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Figure 4. Plots of phase particles on the radial mesh points at the final state in the (/,r)— plane for r} = 0.1 with a = 0.1, 1, and 10, 
where / represents the phase-space density of the particles, and r denotes their radial position. 
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Figure 5. Local maximum pliasc-space density in tlic relaxed state, /m, for r] = 0.1 with a = 0.1, 1, and 10. 
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Figure 6. Fractional mass distribution, dM{j)/dj, against the angular momentum, j, for the initial models. Notice that the mass 
distribution in angular momentum is independent of the initial virial ratio because the initial size of the sphere, Rq, is set to be a 
reciprocal of the virial ratio. 
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Figure 7. Relaxed density profiles for the power-law density spheres of p oc r " with 71=0.5, 1, and 2. The virial ratio and the anisotropic 
parameter are set to be »; = 0.1 and a = 10, respectively. The straight dotted line shows the density profile with p oc r~^'^. 
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